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B physics at CDF - the Beauty of hadron collisions 

D. Tonelli a * 
a Fermilab, 

P.O. Box 500, Batavia, IL, 60510, USA 

The CDF experiment at the Tevatron pp collider established that extensive and detailed exploration of the 
6-quark dynamics is possible in hadron collisions, with results competitive and supplementary to those from e + e~ 
colliders. This provides an unique, rich, and highly rewarding program that is currently reaching full maturity. I 
report a few recent world-leading results on rare decays, CP-violation in B® mixing, and b — ► s penguin decays. 



1. Introduction 

Precise results from the successful B-factories 
experiments disfavor large contributions from non 
standard model (SM) physics in tree-dominated 
bottom meson decays. Agreement with the SM 
within theory uncertainties is also manifest in 
higher-order processes, such as K°-K or B°- 
B flavor mixing. The emerging picture confirms 
the Cabibbo-Kobayashi-Maskawa (CKM) ansatz 
as the leading pattern of flavor dynamics. Non- 
SM contributions, if any, are small corrections, 
or appear well beyond the TeV scale (and the 
LHC reach), or have an unnatural, highly fine- 
tuned flavor structure that escaped all experi- 
mental tests to date. The last chances to avoid 
such a disappointing impasse include the physics 
of bottom-strange mesons, still fairly unexplored, 
along with a few rare B° decays, not fully probed 
at the B-factories because of limited event statis- 
tics. 

The CDF experiment at the Tevatron is cur- 
rently leading the exploration of this physics, ow- 
ing to CP-symmetric initial states in s/s = 1.96 
TeV pp collisions, large event samples collected 
by a well-understood detector, and mature anal- 
ysis techniques. CDF is a multipurpose mag- 
netic spectrometer surrounded by 4ir calorime- 
ters and muon detectors. Most relevant for B 
physics are the tracking, particle-identification 
(PID) and muon detectors, and the trigger sys- 
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tern. Six layers of double-readout silicon mi- 
crostrip detectors between 2.5 and 22 cm from 
the beam, and a single- readout layer at 1.5 cm ra- 
dius, provide precise vertex reconstruction, with 
approximately 15 (70) /jm resolution in the az- 
imuthal (longitudinal) direction. A drift cham- 
ber provides 96 samplings of three-dimensional 
charged-particles trajectories between 40 and 140 
cm radii in < 1, for a transverse momentum 
resolution of <j Pt /p 2 T = 0.1%/(GeV/c). Specific 
ionization measurements in the chamber allow 
1.5(7 separation between charged kaons and pi- 
ons, approximately constant at momenta larger 
than 2 GeV/c. A comparable identification is 
achieved at lower momenta by an array of scin- 
tillator bars at 140 cm radius, which measure the 
time-of flight. Muons with p T > 1.5(2.2) GeV/c 
are detected by planar drift chambers at |ry| < 0.6 
(0.6 < \r)\ < 1.0). 

Low-px dimuon triggers select J/ip, rare B, 
and bottomonia decays. They collected approxi- 
mately 40 millions J/tp decays in 5 fb _1 of data, 
used to reconstruct roughly 6,000 B® — > J/tp(j), 
20,000 B° -> J/V'iC*(892)°, and 52,000 B+ -> 
J/i/;K + decays. A trigger on charged parti- 
cles displaced from the primary vertex collects 
hadronic heavy-flavor decays. It relies on dedi- 
cated custom electronics to reconstruct tracks in 
the silicon with offline-like (48 /xm) impact pa- 
rameter resolution, within 20 fis of the trigger 
latency. This yielded approximately 50 millions 
D° -> K-7T+, 13,000 B° -> D s -(7r+7r-)7r+, and 
12,000 B" -> A'+tt" decays in 5 ftr 1 of data. 
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CDF has currently collected 8 fb 1 of physics- 
quality data. The sample size will reach 10 fb in 
October 2011. Additional 6 fb -1 will be col- 
lected if the proposed three-year extension will 
be funded. 

In the following I report some recent, world- 
leading results, selected among those more sensi- 
tive to the presence of non-SM particles or cou- 
plings. Branching fractions indicate CP-averages, 
K*° is shorthand for the if* (892)° meson, and 
charge-conjugate decays are implied everywhere. 

2. Polarization in B® —> (fxfi decays 

The B® — > (jxj) decay proceeds through a 
penguin-dominated b — > sss transition and was 
first detected by CDF a few years back. We 
recently improved the measurement of branch- 
ing fraction using 295 ± 20 events reconstructed 
in 2.9 fb -1 of data collected by the displaced- 
track trigger. We obtain B(B° S 00) = [2.40 ± 
0.21 (stat.) ± 0.27 (syst.) ± 0.82(6)] x 10~ 5 pQ. 
The last uncertainty is due to the uncertainty on 
B(B° — > J/ip(j>), used as a reference. 

Three polarization amplitudes enrich the phe- 
nomenology of the -B° — >■ 00 decay. They cor- 
respond to the allowed values of orbital angu- 
lar momentum (£ = 0, 1, or 2) in the decay of 
a pseudo-scalar into two vector particles (B — > 
VV). A determination of these amplitudes may 
contribute useful insight into the puzzling pic- 
ture of B VV polarizations, where experi- 
mental data in £?°M — >• <j)K*^ decays disfavor 
first-order predictions of small transverse compo- 
nent (/t = 0(my /rn 2 B ) « 4%, where m are the 
masses), possibly suggesting contributions from 
non-SM amplitudes [2J. 

We report the first measurement of B® — > 00 
polarization, using the sample employed for the 
branching fraction measurement [3]. We fit the 
mass of the four kaons and their angular distri- 
butions (in helicity basis). The angular accep- 
tance, extracted from simulation, is validated by 
measuring the polarization of B® — > J/ip<fi decays 
to be consistent with results from an indepen- 
dent sample [4]. The small bias due to decay- 
length requirements in the trigger is included in 
the systematic uncertainties. In analogy with 



measurements of similar b — > s penguin decays, 
the results are at odds with naive theory predic- 
tions. We measure \A \ 2 = 0.348±0.041 (stat.)± 
0.021 (syst.), |A||| 2 = 0.287 ± 0.043 (stat.) ± 
0.011 (syst.), \A ± \ 2 = 0.365 ± 0.044 (stat.) ± 
0.027 (syst.) for the polarization amplitudes, and 
cos((5||) = — 0.91±^| (stat.) ±0.09 (syst.) for the 
phase difference between the A|| and Aq am- 
plitudes. Further extensions of this analysis to 
larger samples will provide information on the 
decay-width difference Ar in B® mesons. 

3. Rare B — > /J + ^t~ and B — > hfj, + n~ decays 

Decays mediated by flavor changing neutral 
currents, such as B® s ^ —} fx + or B — > h[i + 
are highly suppressed in the SM because they 
occur only through higher order loop diagrams. 
Their phenomenology provide enhanced sensitiv- 
ity to a broad class of non-SM contributions. 

The B®^ — > rate is proportional to the 

CKM matrix element |Vtd| 2 (|Vt s | 2 ), and is further 
suppressed by helicity factors. The SM expecta- 
tions for these branching fractions are O(10~ 9 ), 
ten times smaller than the current experimental 
sensitivity. An observation of these decays at the 
Tevatron would unambiguously indicate physics 
beyond the SM. Or, even improved exclusion- 
limits strongly constrain the available space of 
parameters of several SUSY models. 

The latest CDF search for B9 % — > fJ, + /i~ decays 
uses 3.7 fb -1 of data collected by the dimuon 
trigger with Pt(m) > 2 GeV/c [5]. A loose pres- 
election based on opposite-charge dimuon trans- 
verse momentum and muon identification criteria 
(quality of matching with track, energy deposit in 
the calorimeter, specific ionization energy loss), 
rejects combinatoric background and charmless 
B decays. Further rejection is achieved by se- 
lecting on the decay-length significance (\/a\) 
against prompt background and on the isolation 
of the B?s meson candidate, to exploit the harder 
fragmentation of &-mesons with respect to light- 
quark background. In addition, we require the 
candidate to point back to the primary vertex 
to further reduce combinatoric background and 
partially reconstructed 6-hadron decays. The re- 
sulting event sample contains about 55,000 can- 
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didates, mostly coming from combinatoric back- 
ground. To further enhance purity we use an ar- 
tificial neural network classifier (NN) that com- 
bines information from the above observables into 
a single scalar discriminating quantity. Signal dis- 
tributions are modeled from detailed simulation; 
backgrounds from mass-sidebands in data. The 
signal rate is normalized to 20,000 B + — > J/'0(— > 
H + H~)K + reconstructed in the same sample. The 
ratio of trigger acceptances between signal and 
normalization modes (~ 25%) is derived from 
simulation, the relative trigger efficiencies (~ 1) 
are extracted from unbiased data and the rel- 
ative offline-selection efficiency (~ 80%) is de- 
termined from simulation and data. The ex- 
pected average number of background events in 
the search region is obtained by extrapolating 
events from the mass-sidebands. The validity of 
this extrapolation is checked by comparing pre- 
dicted and observed background yields in sev- 
eral independent control samples including like- 
sign dimuons, opposite-sign dimuons with neg- 
ative decay-length, and opposite-sign dimuons 
with one muon failing the muon-quality require- 
ments. Small contributions of punch-through 
hadrons from B? ^ — > h + h~ decays are included 
in the estimate of total background. We optimize 
the limit by combining three independent ranges 
for the NN discriminator (with efficiencies rang- 
ing from 44 to 12%) to obtain the a priori best ex- 
pected 90% C.L. upper limit on B(B9, ^ + A*~) 
(see fig. Q}. The resulting 90 (95)% CL upper- 
limits are 

B{B° -> ft+fT) < 3.6(4.3) x 10~ 8 , (1) 
B(B° -> fi+fi-) < 6.0(7.6) x 10~ 9 . (2) 

These results are the most stringent currently 
available and reduce significantly the allowed pa- 
rameter space for a broad range of SUSY models. 
An update of this analysis with approximately 
doubled sample size will further improve them 
soon. 

We also updated the analysis of B° -> K*°(-> 
K + 7r-)fi+fj,-, B+ -> K+n+V, and B° -> 
K+K-)fj,+fi- decays to 4.4 fb^of data [6]. 
These are suppressed in the SM {B ~ 10~ 6 ), 
with amplitudes dominated by penguin and box 
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Figure 1. Distribution of mass in three 

independent ranges of the NN classifier. 



b — >• s transitions. Despite the presence of final- 
state hadrons, accurate predictions greatly sensi- 
tive to non-SM contributions are possible for rel- 
ative quantities based on angular-distributions of 
final state particles. The B + and B° decays have 
been previously studied, while the B® channel has 
not been observed yet. We use a large sample 
of high-purity dimuon candidates combined in a 
common vertex with one or two charged particles, 
whose masses are assigned as appropriate for each 
final state. In the B° decay, the mass hypothe- 
sis yielding the Kit mass closer to the known K*° 
mass is chosen, resulting in the proper assignment 
about 92% of the times. A NN that uses informa- 
tion from kinematics and particle-identification 
greatly improves signal-to-background discrimi- 
nation. It is trained on simulated signals and 
mass-sideband data. The simulation is tuned to 
reproduce accurately the data using the corre- 
sponding resonant channels (B — > J/iph). The 
NN is optimized for maximum expected statisti- 
cal resolution on the branching ratio and asym- 
metry measurements. 

Prominent signals of 120 ± 16 B + K + n + fi~ 
and 101 ± 12 B° -> K*°fx+fi- events (fig. [3 left) 
are observed. The absolute branching fractions, 
measured using the resonant decays as a refer- 
ence, are [0.38 ±0.05 (stat.)± 0.03 {syst.)} x 10~ 6 
and [1.06 ±0.14 (siaf.)±0.09 (syst.)} x 10~ 6 , re- 
spectively, consistent and competitive with pre- 
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Figure 2. Distribution of K + n~ /i + [i~ mass (left) 

and K + K~ (i + fi~ (right). Figure 3. Forward-backward asymmetry in B° 

K*° ji + n~ as a function of dimuon mass. 



4. Measurement of the B a s mixing phase 



vious determinations [2]. In addition, 27 ± 6 
Bg events are reconstructed, corre- 

sponding to the first observation of this decay, 
with a significance in excess of 6a. The branch- 
ing ratio, (1.44±0.33 (sto<.)±0.46 {syst. )) x 1CT 5 , 
is consistent with theory predictions, and corre- 
sponds to the rarest B® decay ever observed to 
date. We use the B + and B° signals for the first 
measurement in hadron collisions of branching ra- 
tios, muon forward-backward asymmetry (Afb), 
and K*° longitudinal polarization, as a function 
of the dimuon mass. The asymmetry is greatly 
sensitive to non-SM particles and is determined 
from a fit to the cos(0 M ) distributions, 9^ being 
the helicity angle between the fi + and the 

opposite of the B (B) direction in the dimuon 
rest frame. Angular acceptances are determined 
from simulation. Figure [3] shows the asymme- 
try as a function of dimuon mass. Integrated 
in the 1-6 GeV 2 /c 4 range, where theory predic- 
tions are most reliable, it equals Afb(1 < q 2 < 
6) = 0.43;£ ^ (stat.) ± 0.06 (syst.), consistent 
with Belle and Babar determinations [2]. CDF 
plans to achieve world-leading results by summer 
2011 with the anticipated 2-3 factors increase in 
statistics, due to additional data, triggers, and 
reconstructed final states. 



Non-SM contributions have not yet been ex- 
cluded in B° s -B s mixing. Their magnitude is 
constrained to be small by the precise determi- 
nation of the frequency [7]. However, knowledge 
of only the frequency leaves possible non-SM con- 
tributions to the (CP-violating) mixing phase un- 
constrained. The time evolution of flavor-tagged 
B® — > J/ip4> decays allow a determination of 
this phase largely free from theory uncertain- 
ties. These decays probe the phase-difference be- 
tween the mixing and the b — > ccs quark-level 
transition, /3 s J/lW = /3f M + /3 S NP , which equals 
pSM = &rg(-V ts V* b /V cs V* b ) « 0.02 in the SM 
and is extremely sensitive to non-SM physics in 
the mixing. A non-SM contribution (/3^ p ) would 
also enter fyl 1 ^ = 0f M - 2/3f p , which is the 
phase difference between mixing and decay into 
final states common to B® and B s , and is also 
tiny in the SM: 0f M = arg(-Mi 2 /T 12 ) ~ 0.004. 
Because the SM values for $ J J^ and § J J^"^ 
cannot be resolved with the precision of cur- 
rent experiments, the following approximation is 
used: S J/ ^ a -2/3f p » -2/3 s J/ ^, which holds 
in case of sizable non-SM contributions. Note 
that the phase (frs also modifies the decay- 
width difference between light and heavy states, 
Ar = T L - T H = 2\T 12 \cos((j) J s /i " t '), which en- 
ters in the B® — > J/ip<f> amplitude and equals 
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Ar SM w 2|r 12 | = 0.086 ± 0.025 ps" 1 in the SM 

We updated the measurement of the time- 
evolution of flavor-tagged B® —> J/tJ)(—> 
/i + /x _ )(/>(—>• _ftr + AT _ ) decays to a sample of 
5.2 fb" 1 collected by the dimuon trigger .4 . Im- 
provements over the previous version of the anal- 
ysis include (a) a doubled event sample along 
with a newly optimized selection (b) a fully 
data-driven recalibration of the flavor-tagging al- 
gorithms (c) inclusion of possible non-</> scalar 
K+K~ contributions as B° -> J/ipf (980). 




Mass(J/y 4>) [GeV/C 2 ] Mixing Frequency in ps" 1 



Figure 4. Distribution of J/ipK + K mass after 
the optimized selection. 



A kinematic fit to a common space-point is ap- 
plied to the candidate J/ip and two tracks con- 
sistent with being kaons originated from a 4> me- 
son decay. A NN trained on simulated data (to 
identify signal) and _B° mass sidebands (for back- 
ground) exploits kinematic and PID information 
for an unbiased optimization of the selection: we 
minimize the average expected statistical uncer- 
tainty on the mixing phase extracted from en- 
sembles of statistical trials generated with dif- 
ferent input values. Discriminating observables 
used include kaon-likelihood, from the combina- 
tion of dE/dx and TOF information; transverse 
momenta of the B° s and <fi mesons; the K + K~ 
mass; and the quality of the vertex fit. The fi- 
nal sample contains approximately 6500 signal 
events over a comparable background (fig.|H left). 
The angular distributions of final state particles 
provide a statistical determination of the CP- 
composition of the signal thus improving sensi- 



tivity to the phase. Possible contributions from 
scalar J/ipK + K~ or J /tpf q(980) final states are 
included in the angular distributions. The pro- 
duction flavor is inferred using two classes of al- 
gorithms. Opposite-side tags (OST) exploit fla- 
vor conservation in strong interaction. Being b- 
quarks predominantly produced in pairs with b— 
quarks, the production flavor is inferred from the 
charge of decay products (e, /i, or jet) of the b- 
hadron emitted in the opposite hemisphere of the 
signal B°. Same-side tags (SST) exploit flavor 
conservation in the fragmentation process that 
leads to the signal B® meson. Because a ss pair 
will be required to form a meson from the bottom 
quark, the charge of a kaon kinematically close to 
the bottom-strange meson is correlated to its fla- 
vor. The performance of flavor taggers is modeled 
as function of many event properties using simu- 
lations; however overall scale factors between the 
simulation and reality are allowed and extracted 
from data. The uncertainty in their determina- 
tion contributes to the final systematic uncer- 
tainty of the measurement. The tagging power, 
eD 2 , is the product of an efficiency e, the fraction 
of candidates with a flavor tag, and the square 
of the dilution D = 1 — 2w, where w is the mis- 
tag probability. We calibrate independently b and 
b OST tags using 52,000 B+ -> J/tpK+ decays. 
Predicted and measured dilutions agree with scale 
factors close to unity. The observed tagging effi- 
ciency is (94.2 ± 0.4)% and the average predicted 
dilution on signal is y/ (D 2 ) = 0.110 ± 0.002, with 
1.03 ± 0.06 scale factor. The SST calibration is 
obtained by repeating the full mixing analysis 
on 13,000 B° s -> D-ir+(ir+Tr-) \7[. The result- 
ing frequency, Am s = 17.79 ± 0.07 (stat.) ps _1 , 
is fully consistent with published values (fig. |U 
right). The scale factor, A = 0.94±0.15 (stat.)± 
0.13 (syst.) is given by the size of the amplitude 
at the mixing frequency. This corresponds to a 
tagging efficiency of (52.2 ± 0.7)% and an aver- 
age dilution on signal of ^(D 2 ) = 0.275 ± 0.003. 
Multiple tags, if any, are combined as indepen- 
dent for a total tagging power eD 2 rj 4.5%. The 
proper time of the decay and its resolution are 
known on a per-candidate basis with an average 
resolution of approximately 90 fs _1 . Information 
on Bg candidate mass and its uncertainty, angles 
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between final state particles' trajectories (to ex- 
tract the CP-composition), production flavor, and 
decay length and its resolution are used as observ- 
ables in a multivariate unbinned likelihood fit of 
the time evolution that accounts for direct de- 
cay amplitude, mixing followed by the decay, and 
their interference. Direct CP-violation is expected 
small and neglected. The fit determines the phase 
pj/ijjtf) ^ ^ e decay- width difference Ar, and many 
other "nuisance" parameters including the mean 
B° lifetime (2/(T L +T H )), the magnitudes of lin- 
ear polarization amplitudes, the CP-conserving 
phases (5\\ = arg( J 4|| Aq), S± = arg(A^Ao)), and 
others. The acceptance of the detector is calcu- 
lated from simulation and found to be consistent 
with angular distributions of random combina- 
tions of four tracks in data; the fit model was 
validated by measuring lifetime and polarization 
amplitudes in B° — > J/?pK* decays to be consis- 
tent with i3-factories measurements [9]. 

For enhanced precision, all parameters but the 
mixing phase are determined in a fit with phase 
fixed to zero: ct(B°) = 458.6 ± 7.6 (stat.) ± 
3.6 (syst.) /im, Ar = 0.075 ± 0.035 (stat.) ± 
0.010 (syst.) ps" 1 , |A||| 2 = 0.231±0.014 (stat.) ± 
0.015 (syst.), \A \ 2 = 0.524 ± 0.013 (stat.) ± 
0.015 (syst. ),and 5± = 2.95 ± 0.64 (stat.) ± 
0.07 (syst.) rad. These results represent the cur- 
rent best measurements of these quantities from 
a single experiment. 

With phase floating, fits on simulated samples 
show biased, non-Gaussian distributions of esti- 
mates and multiple maxima, because of known 
likelihood symmetries. We use a frequentist con- 
fidence region in the (/3 s J/ ^,Ar) plane using 
a profile-likelihood ratio ordering [10] . which is 
close to optimal for limiting the impact of sys- 
tematic uncertainties. These are included by ran- 
domly sampling a limited number of points in the 
space of all nuisance parameters. A specific value 
(ft W ,AT) is excluded only if it can be excluded 
for any assumed value of the nuisance parameters 
within 5cr of their estimate on data. 

The resulting allowed region is greatly reduced 
with respect to the previous measurement and is 
fairly consistent (0.8cr) with the SM: the range 
[0.02, 0.52]U[1.08, 1.55] contains the /3 S J/,W phase 
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Figure 5. Confidence region in the (/3 s J/w ,Ar) 
plane. 



at the 68% CL and the range [-tt/2, -1.44] U[-0.13, 

0. 681.[0.89, tt/2] at the 95% CL. 

With the full Run II data sample, we expect 
to observe a non-SM phase or exclude it for any 
value of pi larger than approximately 0.4 rad. 
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